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Abstract 


The physical properties of crystalline compounds are highly essential for predicting their 
applications in various nanodevices. Density functional theory is very useful in that regard. 
Herein, some physical properties such as structural stability, and electronic, optical and 
mechanical properties of two-phase bulk WSez have been comparatively investigated by first 
principle calculations. Band structure calculations elucidate that the materials are indirect band 
gap semiconductors having band gap energy in the range 1-1.5 eV. It was observed that the 
hexagonal crystal has higher absorption coefficients compared to the trigonal phase. Analysis 
of elastic constants of the materials carried out using the Born criteria method indicates that all 
systems are mechanically stable at room temperature. Our findings are in close agreement with 
the previous experimental and theoretical results, this further suggests that the materials are 
promising in optoelectronics applications. 
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I. INTRODUCTION 


becca Diselenide (WSe2) is an inorganic crystalline 
compound belonging to the family of transition metal 
dichalcogenide (TMDC)[1] represented with general formula 
MX2, where M stands for transition metal and X is 
chalcogenides (S, Se and Te) [2]. WSe2 is chemically stable 
[3], p-type in nature and has an indirect band gap [4, 5]. It 
takes hexagonal and trigonal phases. Recently research on 
TMDC materials has shown that dimensionality plays a very 
important role in determining the basic physical properties of 
these compounds [6]. In this regard, a monolayer or two- 
dimension, bilayer and a few layers of TMDCs have been 
successfully synthesized [7]. The bulk and few layers of WSe2 
have found applications in many devices and electronic 
components, such as solar cells, rechargeable lithium 
batteries, photodetectors, transistors and light emitting diode 
(LED) [8]. It is also used as a lubricant. Among TMDCs WSe2 
has been a very attractive crystalline compound in recent years 
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due to its precious physical properties such as good electron 
mobility, structural stability, and band gap in the visible light 
spectrum (400 to 700 nm) [9]. Reference [10] reported a GW 
and PBE method to examine the band structure of bulk WSez. 
More recently, [6] investigated the electronic and optical 
properties of hexagonal bulk WSe2 (p63/mmc) using DFT, 
utilizing PBE-GGA. They confirmed the semiconducting 
nature of this compound. Though the physical properties of 
the bulk hexagonal phase of WSe2 are well studied, the 
trigonal phase has received far less attention and its properties 
are not completely explored. 

In this paper, we used first-principles density functional 
theory (DFT) to comparatively investigate the structural, 
electronic, elastic and optical properties of the hexagonal with 
space group (p63/mmc), group number 164 and trigonal-phase 
having space group (p3m1) and group number 187 bulk 
WSe>. 
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II]. COMPUTATIONAL DETAILS 


In this study, we used the projector augmented plane wave 
(PAW) numerical method to solve the Kohn-Sham equation 
within the framework of density functional theory (DFT) 
implemented in Quantum ESPRESSO code [11] to investigate 
the hexagonal and trigonal phases of bulk tungsten Diselenide 
(WSez) with 6 and 3 atoms respectively. First principle 
functional theory non-spin polarized calculations [12] were 
carried out in the generalized gradient approximation for 
crystals as given by Pardew, Burke and Ernzerhof (GGA- 
PBE) [13]. For integrals smearing has been employed and 
precisely the Maxfesse-Paxton smearing method [14], for both 
the two crystalline compounds. Monkhost-Pack Scheme has 
been adopted for the Brillouin zone integration [15] with 
12x12x3 k-points grids for WSe:, in the electronic properties 
calculation. The WSez structures were optimized using 
Broyden-Fletcher-Goldfarb-Shannon (BFGS) algorithms[16]. 
The atomic positions were fully relaxed until the forces acting 
on the atoms were below 0.001 eV/A. 


III. RESULTS AND DISCUSSION 


A. Structural properties 


The equilibrium structural ground state parameters of the 
two-phase bulk WSe2: compounds, such as the lattice 
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constants, the structural optimization in the quantum espresso 
code has been carried out using Perdew Burke Ernzerhof 
generalized gradient approximation (PBE-GGA) by 
minimizing the total energy with respect to the volume of the 
unit cell. The geometric structures of the compounds are 
hexagonal and trigonal as shown in Fig. la and b, Fig. Ic 
shows the Brillouin zone used for k-path in band structure 
calculations for the considered crystal compounds. Fig. 2 
depicts the energy as a function of the volume of the 
hexagonal and trigonal phases of the considered compounds. 
The calculated total energies are fitted to Murnaghan's 
equation of state to obtain the ground state equilibrium lattice 
constant. The total ground state energies for the hexagonal and 
trigonal structures are -106.207316 and -53.103372Ry per 
supercell which confirm the structural stability of the 
considered systems. The calculated and experimental lattice 
constants, bulk modulus and its pressure derivative are listed 
in Table I. It is seen that the lattice constants obtained with 
PBE-GGA agree with other similar calculations. The 
calculated bulk modulus in Hexagonal WSez, 59.7 GPa, is 
15% smaller than that of trigonal. 


(©) 


Fig. 1 Primitive cells of Bulk (a) Trigonal, (b) Hexagonal WSez structures and (c) Brillouin zone used for k-path in band 
structure plots. 


-106.185 4 
-106.190 4 


-106.195 4 


QO 


1550 


Energy (Ry) 


-106.200 4 


| 


1450 1500 


-106.205 


-106.210 


T ¥ T T T 
1600 1650 1700 1750 1800 


Volume (a.u‘3) 


Energy (Ry) 


T T T 
1000 1050 1100 1150 1200 


Volume (a.u*3) 


Fig. 2 Energy volume curves of Bulk (a) Hexagonal and (b) Trigonal WSez structures. 
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Table I. Calculated equilibrium lattice constants for Bulk 
Hexagonal and Trigonal WSez structures by using PBE-GGA. 


Ref. a(A) c (A) B(GPa) _B Eo (Ry) 
Hexagonal WSe2 3.356 12.743 23.15 2.492 -106.20731 
*Hexagonal WSe, 3.216 [10] 13.170 [10] 

3.628 [17] 12.980 [17] 
Trigonal WSe, 3.543 12.876 17.10 2.878 -53.10337 
*Trigonal WSe, --- — 
*PBE GGA calculations 


B. Electronic Properties 


We use the optimized structures of the considered materials 
to carry out the electronic properties’ calculations. In 
electronic band structure calculation, we used PBE-GGA. The 
band structures of the two compounds were calculated along 
high-symmetry directions of the Brillouin zone [-M-K-I° and 
the fermi energy level is set at 0 eV. The calculated band 
structures for the considered crystals are shown in Fig. 3. Both 
phases of bulk WSez are indirect band gap semiconductors. 
The calculated band gap energies are 1.4 and 1.5 eV for 
hexagonal and trigonal respectively, these values together 
with experimental values are listed in Table II. These results 
show that DFT has overestimated the band gaps. This is due 


k-point mesh 
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to the inherent drawback of GGA/LDA functional courses by 
self-interaction errors [18, 19]. Fig. 4 illustrated the density of 
state plots for the considered compounds, they all confirmed 
the band gap energies of the materials as observed in band 
structure plots. 

Electron acquires masses, different from its rest mass (mo) 
when accelerated in a magnetic or electric field. This mass is 
known as effective mass (m*) and has a negative value in the 
valence band and a positive value in the conduction band. By 
using the curve fitting technique, the effective mass of the 
electron in conduction band minimum (CBM) and hole in 
valence band maximum (VBM) were obtained from the E-K 
plot (Band structure) of the compounds by applying (1). 

dk? 


m* = hh? [se 
dE 


(1) 


Where h = h/2z, h is planks constant, k is wavevector and 
E is energy. 

Table III shows the calculated effective masses of electrons 
and holes for the two phases in the unit of electron rest mass 
(mo). Reference [20] has reported small effective masses of 
electron holes for perovskite/HfS2 and attributed it to the 
efficient separation of photon-generated excitons. This is what 
makes them attractive for photovoltaic applications. 
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Fig. 3 Band structure of Bulk (a) Hexagonal and (b) Trigonal WSe> crystalline compounds. 
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Fig. 4 Density of state of Bulk (a) Hexagonal and (b) Trigonal WSez crystalline compounds. 


Table II. Calculate band gap for Hexagonal and Trigonal bulk 
WSe Structures. 


Compound Band Gap (eV) 
This work DFT (PBE-GGA) Experimental 
Hexagonal WSe, 14 1.00 [6], 1.33 [10] 1.75 [10] 


Trigonal WSe, 1s 1.10 [17] 


Table HI. Calculated effective mass of electrons and holes for 
Hexagonal and Trigonal bulk WSez Structures. 
Compound Effective mass (GGA-PBE) 
m*,(Kg) m*,(Kg) 
Hexagonal WSe> 0.104m, 0.137m, 
Trigonal WSe 0.099m, 0.133m, 


*Experimental 


C. Optical Properties 

The crystal's response to incident radiation determined its 
application for optoelectronic devices and nanoelectronic 
components [21]. This prompted the study of the optical 
properties of crystalline materials, such as absorption 
coefficient, optical conductivity, refractive index, reflectivity 
etc. These parameters can be determined from the dielectric 
tensor €(@) given by [22]. 


€(W) = €,(@) + ie, (@) (2) 


The imaginary part €2(w)) can be calculated from the 
Kubo-Greenwood equation expressed as (3)[23], 


2me? 


€2(W) = Fo Lewcl(Pelu.rl Pel? OER — Ee — E) (3) 
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While the real part €;(w) of the dielectric tensor can be 
obtain from the imaginary component using the Kramer- 
Kronig equation given by, 


€,(w) = 1+ @) i dey 229 (4) 


w2—q2 


a(o) = V20| fa)? + f2@F-a@) 6) 


o(w) = (6) 
n= v(HeOE Heo sae) (7) 


Where w is the frequency of phonon, © is the unit cell volume, 
u is the unit vector, e is the electron charge, Y¢ and Yf are the 
wavefunctions for conduction and valence band electrons 
respectively at k, a(w) is the absorption coefficient, o(w) 
optical conductivity and n is the refractive index of the 
material. 

We investigated the optical properties in the incident photon 
energy range of zero to 20 eV. Fig. 5a and b show the plots of 
the real part of the dielectric tensor €)(@) against energy and 
imaginary part €2(@) against energy for the hexagonal system 
respectively. The graphs show that there is no significant 
difference along the x and y-directions for the values of both 
components. This is also observed with trigonal WSe2 
structure Fig. 5c and d, this could be attributed to their lattice 
constants, where a is equal to b but not equal to lattice 
parameter c. The values of the static real part of the tensor 
€,(0) in the x and y-direction are 5.93 while in the z-direction 
is 2.40 for the hexagonal compound, as recorded in Table IV. 
The trigonal phase €1(0) in the x and y-direction has a value of 
4.19, and 2.35 in the z-direction. This confirmed that both 
structures are anisotropic. 
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Table IV. Calculated static dielectric €1(0) for Hexagonal and 
Trigonal WSez bulk structures in three polarization directions. 


Compound €:(0) 

Elx Ely Elz 
Hexagonal WSe. 5.92954670 5.92954670 2.40032769 
Trigonal WSe> 4.19333697 4.19333697 2.34844347 


As we compared to the hexagonal, the dielectric behaviour 
for trigonal shown in Fig. 5 is slightly the same. However, the 
values of the dielectric tensor for hexagonal are much larger 
than trigonal. The zero-frequency limit ¢1(0) shown in the Fig. 
5a gives the static dielectric tensor. The calculated static 
dielectric functions of both systems for different polarization 
directions are given in Table IV, the calculated dielectric 
tensors of the hexagonal are a little bit higher than trigonal. 
Fig. 6a shows the variation of the refractive index and 6b 
absorption coefficient in the energy range from zero to 20 eV 
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for the hexagonal compound. Whereas Fig. 6c and 6d show 
the same for the trigonal phase. In the absorption plot in Fig. 
6b and 6d the first maxima occur at 1.48 eV and 1.53 eV 
corresponding with band gap energies of the materials earlier 
stated in Table II. 


Table V. Calculated static refractive index n (0) for Hexagonal 
and Trigonal WSez bulk structures in three polarization 
directions. 


Compound n(0) 

Elx Ely Elz 
Hexagonal WSe, 2.43785851 2.43785851 1.55387872 
Trigonal WSe. —2.04971319 2.04971319 1.51912046 


Fig. 6c and d depict the plots of refractive index with photon 
energy. From these graphs, we extract the static refractive 
indices (n (0)) of the two-phase WSez as seen in Table V. 


Photon energy (eV) 


(d) 


Fig. 5 Dielectric tensors of Bulk (a and b) Hexagonal and (c and d) Trigonal WSe crystalline compounds. 
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Fig. 6 Absorption (0 (@)) and refractive index (n(@)) plots of Bulk (a and b) Hexagonal and (c and d) Trigonal WSez2 
crystalline compounds. 


D. Mechanical Properties 


Elastic constant and elastic moduli provide the ground for 
understanding the mechanical properties of crystalline solids 
such as brittleness [12], stability, hardness, ductility, and other 
mechanical characteristics. The calculated elastic constant for 
the considered bulk structures within the framework of PBE 
GGA functional is presented in Table VI. These results are 
compared with the theoretical and experimental findings. The 
mechanical stability of bulk materials under isotropic pressure 
can be analyzed using the popular Born stability criteria [9] 
given as, 


C11 > 0,C33 > 0,C,, > 0 (8) 
(Cy — Cy2) > 0 (9) 
(Cy, + Cy2) > Ca3 — 2C7 (10) 
The independent elastic constant is given by, 

Coe = (Cia — C12)? (11) 
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All the studied materials satisfy the Born mechanical stability 
limitations as illustrated by (8)-(11). Furthermore, the table 
demonstrates that the elastic constants for all the materials are 
all positive indicating that they are all stable [9] [24]. 


Table VI. Elastic Constant of Bulk Hexagonal and Trigonal 
WSe within GGA-PBE Functional. 


Elastic constant (Kbar) GGA-PBE 
Hexagonal Trigonal 
WSe2 WSe2 

Ci 801.996 587.88043 
Cr 172.489 117.82029 
Cis 13.858 23.48457 
C33 13.738 55.07496 
C44 21.988 10.92863 
Cee 396,279 220,956 
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1) Hardness and Elastic moduli 

The mechanical properties of a material can be determined 
from its calculated elastic constants. The machinability of a 
material can be calculated by dividing its bulk modulus by its 
elastic constant C44. It was established that the higher the value 


of the ratio (=) the better the machinability [9]. Herein, we 
44, 


have calculated the bulk modulus (B), the shear modulus (G) 
and the Young’s modulus (Y) for the individual C;; for all the 
considered materials. According to Pugh's rule, the ratio (G/B) 
can be used to determine the brittleness or ductility of the 
material compound. This rule states that if < > 0.57 then the 
material is ductile else it behaves as brittle. In this work, the 
calculated values of _ indicate that all the materials have 


shown a ductile nature. Reference [25] reported that material 
ductility can be determined using Poisson's ratio. They 
elucidated that Poisson's ratio greater than 0.25 means the 
material is ductile, and anything less than 0.25 indicates a 
brittle material. 

For hexagonal structures, the bulk (By) and shear (Gy) 
moduli can be obtained from (12) and (13) using the elastic 
constant, 


By = =[(2Cy1 + Cy2) + C33 + Cys] (12) 
Gy = =[C1 + Cy2 — 2033 +4Cy3 +12Cy, +Co5] (13) 
Synthesized systems are more often studied in the form of 
ceramic polycrystalline, hence, it is important to calculate the 
elastic parameters of the investigated materials systems in 
their polycrystalline state. To that effect, we use the Vogt- 
Reuss-Hill approximation, according to which two other 
essential approximations are utilized, i.e., the Voigt and Reuss 
approximations. The expressions for the shear and bulk 
moduli are given by the Reuss model as 
_ (C11 +C13)C33-2(C13)? 


Br Cy +Cy2+2C 33-413 ne) 
6, = 5 ((C1a+12)C33-2(C13)?)Ca4Coe (15) 
R 2 [3ByCaa Coot {(C11+C12)C33-2(C13)*}(Caa+ Coo) 


The effective modulus of anisotropic polycrystalline 
crystals is calculated using bulk modulus and the shear 
modulus, as approximated by Hill [26] is given by, 


G ==(G, + G,) (16) 
B == (B, + By) (17) 
The average and the actual effective moduli above can be 


used to obtain bulk and shear moduli as approximated by 

Voigt-Reuss-Hill given by 
_ 9BG 18 
~ 3B+G (18) 
Now, the Poisson’s ratio(v) for polycrystalline material can 


be given by the expression, 
3B-2G 


~~ 2(3B+G) 
2) Crystal Anisotropic Index 
The orientation dependency of the elastic moduli is called 
crystal elastic anisotropy [27, 28]. Anisotropy of a material 
can be calculated using elastic constant C44 as shown in (20) 


(19) 
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[9]. The anisotropy index has been assigned to crystals for its 
importance in the design of crystals in science and 
engineering. Reference [9] reported that Ag, = 1.0 shows an 
isotropic crystal, while Ag, > 1.0 or Agn < 1.0 indicates 
anisotropy. In this work, the results presented in Table III 
suggest that all our crystals are anisotropic by nature. 


4C44 
Aq = > 
CagtCqgtCyg 


(20) 


IV. CONCLUSION 


This paper reports a computational description and 
comparison of structural, electronic, optical and 
mechanical characteristics of hexagonal and trigonal 
phase WSez bulk systems using DFT calculations. Our 
results elucidate shows that the two Polymorphs of the 
crystal WSez are stable, they are both semiconductors 
with an indirect band gap between 1-1.5 eV, which 
agrees well with previous theoretical and experimental 
findings. The results show that the bandgap of the 
trigonal is larger than hexagonal. The complex dielectric 
tensor, absorption spectrum and refractive index were 
calculated and discussed. Our findings further reveal that 
the hexagonal's real and imaginary parts are larger than 
the real and the imaginary parts of the dielectric function 
of trigonal. The hexagonal has a higher absorption 
coefficient and refractive index in all three polarization 
directions. The results further suggested that the two 
materials are anisotropic and brittle. 
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